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Abstract 

Astrophysical constraints of new physics are often limited to weakly interacting 
light particles, such as axions, the Kaluza-Klein (KK) gravitons from the ADD model, 
sterile neutrinos and unparticles. We discuss the possibility for an astrophysical 
scenario to (dis)confirm new physics for heavy particles beyond TeV energy scale. In 
our scenario, the KK protons (the KK excited quarks/gluons within protons) within 
the framework of universal extra dimensions (UEDs), are produced by high energy 
p + p collisions in Fermi accelerated environments, with protonic isotropic spectrum 
dN/dE oc E^'^ up to at least 10^^ eV. Thus, because they are also electrically 
charged, they should be re-accelerated by mechanism similar to normal protons. The 
KK states (no matter whether they have already decayed to the lightest KK particle 
or not) should contaminate 10~^ to 10~^ of cosmic-ray events for some fixed energy 
E (within some suitable assumptions). Hence, if we have techniques to identify them 
from air shower data, we can constrain UEDs scenario. Our method is an "existence 
proof" that we can constrain new physics beyond TeV scale or much higher by 
classical astrophysical scenarios, which can also be generalized to supersymmetric 
models, the bulk Standard Model fields within the RS model, and the endlessly 
emerging new models. Moreover, it can exploit domains which have no possibility to 
be studied in terrestrial experiments. 
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1 Introduction 



Brane-world scenarios, such as the Arkani-Hamed-Dimopoulos-Dvah (ADD) model [HEIE] 
and Randall- Sundrum (RS) model [H [5], give an alternative framework for solving the 
hierarchy problem. Within the ADD model, the Standard Model (SM) fields are confined 
to a 3-brane ((3 + l)-dimensional spacetime) while gravitons propagate freely in a torus 
compactified bulk space (large extra dimensions), and the gravitational coupling constant 
G = 1/Mpj observed in our (3 + l)-dimensional world is an effective one. So it is natural to 
understand why G is so small. The RS model solves the same problem by a slice of ArfS's 
spacetime. Some superstring-inspired descriptions make these models more attractive. 

Some lineage scenarios of the ADD model also allow the SM particles propagating, 
to some extent, in the bulk spaces. Beside gravitons, such kind of SM particles can also 
be Kaluza-Klein (KK) excited in the extra dimensions, thus give abundant physical phe- 
nomena. A natural extension of the ADD model is let the brane with a finite thickness 
and complex substructures [6]. An example is given in Ref. [7]. In this scenario, quarks 
and leptons are confined to different branes, while the Higgs and SM gauge fields are sand- 
wiched in, hence the possibility for proton decay can be exponentially suppressed. The KK 
excited gauge particles (which may be called branons) can be either baryophobic or lepto- 
phobic, because they feel nontrivial on the brane substructures. Another natural extension 
is asymmetrical compactification, which has two (as minimum, maybe more) separate com- 
pactification scales. The "very large" extra dimensions let only graviton propagate, just as 
what the ADD model says, but the TeV~^ scale "large" extra dimensions, may have the 
SM fields extending [8J. Of course we can on the other hand keep all compact dimensions at 
TeV~ scale rather than make the compactification asymmetric. However, the advantage 
of solving the hierarchy problem in the ADD model is bereft. 

In the case of that the SM fields also extend to some extra dimensions, to obtain 
chiral fermions in the 4-dimensional effective theory, we have only two ways to go: (i) 
to confine fermions to branes only J9], or (ii) to impose bulk fermions orbifold bound- 
ary conditions [lOj. Universal extra dimensions (UEDs) [TT] scenario is an example of 
the second approach. In UEDs scenario, all SM fields can propagate in these "univer- 
sal" extra dimensions, and conservation of momentum in the universal dimensions turns 
to conservation of the KK number in our (3 + l)-dimensional world. For two or more 
universal extra dimensions, the naive KK mode sums diverge when the KK tower number 
iYxK —>■ oo. So, let us consider only one universal extra dimension as minimal univer- 
sal extra dimensions (MUEDs) in this context. In this case, S^/7j2 orbifold compact- 
ification is assumed, and the KK mass eigenvalues have a simple form = n/R, 
where R is the compactification scale. For the reason that only loop diagrams can con- 
tribute electroweak observables by the KK number conservation, the experimental bound 
for UEDs is only M|^^ = > 300 GeV. In the tree level, the mass spectrum of the 
KK excited SM particles has the form MsM,n = \/{M^^yl^M§^, where Msu is the 
zero-mode on-shell mass of the corresponding SM particles. Hence, their masses are 
level- by-level highly degenerated when Msm ^ M^^. However, when radiative cor- 
rections are concerned [121 the mass degeneration is broken, to some extend, as 
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Mg,„ > Mq,„ > Mq,„ > Mw±,n ~ Mzo^„ > Ml,„ > Mi,„ > M^,„ ~ M„, where ^ de- 
notes gluon, Q (L) denotes weak-doublet quark (lepton), q (/) denotes weak-singlet quark 
(lepton), and 7 denotes photon. The KK number conservation breaks down to a KK parity 
that the even and odd KK numbers cannot transform to each other. The correction scale 
depends on some unknown parameters; however, some reasonable choice of parameters 
shows that the largest correction AM^^^ may be as large as 10% [13]. Hence, heavier KK 
excited states should cascade decay (by the KK conserving or even violating interactions) 
to the lightest KK particle (LKP) 71 which is stable [H], by emitting soft SM particles. 
When considering the possibilities of experimental discovery, one always assumes that the 
lifetimes of heavier KK excited states are sufficiently short, thus the states can decay within 
the collider; however, it is not supposed to do so. The total width cannot in fact be cal- 
culated by the failure of reconstructing the Breit-Wigner resonance [15j. Notice that the 
decay rate of an unstable particle dV oc l/m_4 in the phase space formula, the lifetime 
r oc m_4 where m_4 is the mass of the decay particle. If the AMs are smaller for a set 
of parameters different from in [13j, or soft SM cascade processes are suppressed by other 
reasons, the lifetimes of heavier KK excited states should be even longer. Specific calcula- 
tions for whether the not-very- short- lived KK excited states can affect other more mature 
scientific scenarios, such as disturb predictions of Big-bang nucleosynthesis, or distort the 
cosmic microwave background, are need; however, they are beyond our scope of our paper. 
Some naive considerations show that all of them are not very crucial, because we do not 
really need longer-lived KK excited states (although the long-lived ones are also possibili- 
ties we shall consider in the identification section in but some not- very-short-lived KK 
excited states to suffer the time scale of Fermi acceleration (which is maybe ~ s or much 
shorter), which is much shorter than the time scale of the scenarios we mentioned above. 
So we assume that the lifetimes of heavier KK excited states are long enough to suffer the 
astrophysical scenario we draw in this paper. 

Astrophysical constraints of new physics are often limited to weakly interacting light 
particles, such as axions [l6l[I7l[T8l|19l[20l[2ll|23[23l[2l[25ll^ gravitons 
from the ADD model [3l|28l[29l|30l[3ll[32l[33l[3l], sterile neutrinos [351 |36l ISH [38] 
and unparticles [SH]- We want to construct an astrophysical scenario to (dis)confirm new 
physics for heavy particles beyond TeV energy scale. Notice that in Fermi accelerated 
environments, protons in a power law spectrum up to at least lO^'^eV should be produced 
(even if we have already derived an overall Lorentz factor F ~ 300), thus p + p collisions 
up to a tremendously large energy should happen there, which we cannot even imagine in 
terrestrial experiments. However, we have to brain storm to know their happenings. In this 
paper, we construct a scenario which may (dis) confirm UEDs by high energy observation 
of cosmic-rays. This scenario may or may not have opportunities to give stronger bounds 
than colliders, because of the large uncertainties in our estimations, and the technical 
details of lots of synergic scientific domains (which we cannot discuss at length in this 
paper). However, it is at least an "existence proof" for this kind of methodology. It can 
also explore domains which have no possibility to be studied in terrestrial experiments. 
Some similar scenario in Ref. [10], also suggested the production of some kind of strongly 
interacting massive particles mp + p collisions in astrophysical environments; however, our 
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scenario have a lot of advantages than theirs. The advantages rise mainly because (i) one 
of the protons in their scenario stays at rest, but both of the protons in our scenario are 
Fermi accelerated, and (ii) our KK excited states suffer an additional accelerated process. 
Detail comparisons are given in §6) 

In our scenario, the KK protons (with either KK excited quarks or gluons in it) are 
produced hj p + p collisions in Fermi accelerated environments. Both the original Fermi 
mechanism or diffusive shock accelerating model have an isotropic spectrum dN/ dE oc E~'^ 
up to at least 10^^ eV, hence they are okay for our purpose. The KK protons should also 
be accelerated just as normal protons by the same mechanism; however, they should have 
different properties than normal ones. Beside being discovered one by one from air shower 
data directly, they may be accelerated to energies normal protons cannot be accelerated 
to, or they (or their decayed final state) may contaminate significant amount of ultra- 
high-energy cosmic-ray (UHECR) events because of an overall energy shift, both of which 
may make them a discovery. In ^ we calculate the cross section and production rate 
of the KK protons. We show that the production rate may be large enough to make 
meaningful scientific constraints. In [|3|, we discuss the accelerated property of them in 
Fermi accelerated environments, making a comparison with normal protons. We notice 
that the KK states should contaminate 10~^ (for special sources) to 10~^ (for diffuse flux) 
of cosmic-ray events for some fixed energy E (if assuming the optical depth Tpp = 1), which 
are not too small a sample to be discovered by air shower detection. We also discuss the 
propagating properties of them related to soft photon interactions. In §11 we consider the 
probabilities to identify them (or their decayed final state) from other cosmic-ray particles 
from air shower data. If it can be done so, our method can have larger possibilities to give 
smaller parameter space for UEDs than other methods. In ^ we calculate the possible 
constraints of the KK cosmic-ray flux from neutrino detectors; however, the constraints are 
very loose for current scientific equipments to affect our former estimations. We discuss 
our results and draw the possibilities to generalize our method to other new physics models 
in P 

2 Producing of the KK Protons 

The KK particles can be produced by high energy collisions, both in terrestrial experiments 
and astrophysical environments. We may want to consider a list of some not-very-exotic 
processes, such as 7 + 7, e + e, p + p, p + n, e + •y and p + 7. Collisions including 
electrons are very tempting; however, electrons always have limited energy because of 
synchrotron radiation. Ultra-high-energy (UHE) photons with energy 10^^"^'* eV or larger 
may be produced in astrophysical scenarios ^41j; however, we leave behind this possibility 
elsewhere because of page limitations. Collisions including neutrons are very interesting; 
however, uncharged neutrons are always less energetic than protons, thus they make a 
less center-of-mass energy ^/s for p + n interactions than p + p ones. In addition, we lack 
trustworthy parton distribution functions (PDFs) f{x,Q) for high Q neutrons. Hence, we 
will only discuss p + p collisions in this paper. 
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2.1 p + p^ (the KK states) Cross Sections 

The p + p cross sections to produce the KK bosons and fermions in the framework of UEDs 
have aheady been calculated in [l2l HS]. However, their motivation is mainly on whether 
we have possibilities to confirm them on synchrotrons (especially Tevatron I, II and the 
Large Hadron Collider (LHC)), thus they focused their calculations on a center-of-mass 
energy v/i of 1 TeV, 2 TeV and 14 TeV respectively and let Mf ^ ClS db free parameter. As a 
result, their calculations are not suitable for our purpose. So we redo the calculations for a 
set of parameters ^/s. In this paper, we always fix Mf^^ = 1/R = 350 GeV for simplicity. 

In our calculation, we assume the KK states are sufficiently stable (see the discussions 
in We use the amplitude-squared calculated by Macesanu et al. [13]. Their 

expression for total cross section is 



where S" is a statistical factor to memorize the number of identical final states, s, t{u) = 
— (s/2)(l — M^/s)(l =Fcos^^) are the Mandelstam variables, 6 is the angle between two inci- 
dent particles in the source comoving frame. We evaluate the PDFs f{x, Q) by CTEQ6.6m 
(standard MSbar scheme) [HI HS] . The calculation of total cross section a™ is too cum- 
bersome although routinely, because there are too many subprocesses and too many quark 
flavors. Notice that (when Q = 350 GeV) u-quarks dominate in the large x region and 
gluons dominate in the small x region. So we neglect the subprocess for initial states s, c, 
b and anti-quarks completely in the context. We also neglect the contribution of d-quarks 
because their PDFs are always smaller than u-quarks; however, this reduction may not 
be suitable for arbitrary Q. Hence, we will only calculate subprocess including u-quarks 
and gluons. In our calculation, MsM,n = is assumed for simplicity, for the reason 

that both zero-mode on-shell mass and radiative corrections cannot alter the masses of 
magnitude. 

The cross section is calculated in Fig. [H We integrate out the / d9 to neglect the 
transverse momentum distribution, which is irrelevant to astrophysical scenarios. We use 
only the lowest order (LO) expression for quantum chromodynamics (QCD) running cou- 
pling constant as{Q) in the amplitude-squared, which is sufficient for high energy regions 
Q ^ Mzo. We also assume Q = Mn for simplicity. We do not show the uncertainties 
of the PDF sets in that figure; however, some tentative calculations show that error bars 
can be neglectable for a large range of Q, and it will never affect our rough astrophysical 
estimations. As seen from the figure, the contribution of n = 1 state is always more impor- 
tant than higher excited states in the calculable regions, thus we neglect the others. The 
subprocess g + g gn + On dominate the large ^/s region, and q + q ^ Qn dominate 

^However, this assumption may need careful considerations when ^/s is much larger than M„. We 
assume it because a refined calculation is far beyond our scope in this paper. 




(1) 
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the small y/s region. In magnitude, our results are consistent with [13] in the regions not 
adjacent to (2M„)_|_. However, for energy just above the threshold y/s = (2M„)_|., our 
as are largely suppressed rather than mildly decreased in theirs. In addition, it seems 
that their cr^^{g + g — >) increases more softly than ours. These may only be caused by 
some numerical details, so we will not discuss these issues deeply. Hence p + p collisions 
can produce the KK protons with the KK quarks or gluons in it (we will call them by 
a joint name Pkk)- The reason why the subprocess g + g ^ Qn + On ^Iso suitable for 
our purpose to produce the KK excited protons is that in the high energy collisions gluon 
bremsstrahlung is unsuppressed, thus the two gluons will produce two jet events to form 
two protons. Because we have already assumed that the KK number conservation is exact 
in UEDs scenario, a KK excited gluon cannot take off its KK number and cannot leave the 
proton for color confinement. Thus both the KK excited bosons and fermions can form 
the KK excited protons. 




Figure 1: Cross section for the case n = 1, 2, 3 and subprocess g + g ^ gn + 9n^ u + g —>■ 
< + M + M < + <, M + u' < + u'^* and u + u < + < (Eq. (15), (17), 
(18), (22) and (24) in Ref. [13] respectively), as marked in the figure. The total cross 
sections for different n are showed by thick black lines. We have assumed q = u because 
u-quarks dominate the quark PDFs in a large range of energy scale Q. The other process 
g + g ^ (ln + (in be important; however, we neglect it because of a lack of consistency 
amplitude- squared. 

The calculated cross section result has a lot of cumbersome difficulties, mainly in the 
high a/s regions (it has to apply) up to at least 10^^ eV. The difficulties are listed as follows: 
(i) Because the PDFs are only apphcable to pn = 'iM'^/ s > x^i^i = 10~^ in CTEQ6.6 series. 
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we can only calculate cross sections below energy ^Smax = 7 x 10^^ eV d for the lowest 
KK state Mf^^ = 350 GeV, which is insufficient for astrophysical events that happen in 
Fermi accelerated environments, (ii) Multi-KK processes (just like the multi-vr processes 
in soft hadron physics) have to be considered when y/s is much larger than M„. (iii) The 
Lagrangian we use to calculate the cross sections may only be an effective one (thus non- 
renormalizable in the framework of quantum field theory), hence cannot be applied to ^/s 
larger than some cutoff energy scale (which may be not very larger than Mi ) . Problem (i) 
and (ii) are the main topics of §2.2[ Problem (iii) is much more troublesome, thus we can 
only solve it completely after the emergence of the final theory (which should explain what 
is the physical essence of normalization). We assume their capabilities simply because (i) 
some smoothness considerations, and (ii) we do not have other ways to go. Other authors 
used the same strategies as ours, just hke what in Ref. [l6j. Moreover, because of the re- 
acceleration of our KK particles, the energetic KK ones need not to be produced by high 
a/s p + p collisions; hence, the high ^/s behaviors of our cross sections may in fact not very 
crucial. It should be that case at least for the naive estimations we choose in this paper, 
because we can hardly distinguish the difference between the KN and DL parametrizations 
in Fig. m however, we do not know the applicability of this supposition. We leave the more 
quantitative discussions in the follow-up publications. 

2.2 High Energy Behavior 

As we have already mentioned before in §2.11 (i) the maximum energy a numerical total 
cross section cr^f (^Smax) can achieve, depends on the minimum x PDFs /(x, Q) can give. 
For CTEQ6.6 series, the believable maximum energy for the KK ffist-exist state Mf^ 
is only 7 x 10^^ eV, which is insufficient for astrophysical purpose. Moreover, (ii) the 
contributions of multi-KK producting processes lack of estimations in our cross section 
calculations. Naively, one may want to continue cr^{^/s) by fitting any functions putting 
by hand, with our low energy numerical results calculated above in Fig. [H However, if we 
assume that a™ is dominated by some special subprocess (as overpowered g + g ^ gl+gl 
shows in Fig.[T]in our cases), it may globally be inconsistent with the (axiomatic) Froissart- 
Martin bound □ atot < Cln^s for sufficiently large s [17J. Of course, everything will be 
fine, assuming that Froissart-Martin bound is far from being saturated. 

The high energy behavior oi p + p total cross section is questioned in length [18] . In 
this context, continuation of <y^^{y/s) to higher energy will be done by the Kang-Nicolescu 

^We cannot saturate the bound because of our limited computation power. However, increasing the 
calculation load considerably can have little improvement on the result of i?max, thus we do some qualitative 
analyses instead in ^2.2\ 

•^Of course, the Froissart-Martin bound need not to be obeyed for a KK mass tower. Because its 
assumptions include that all masses are equal to the unit of mass when the energy variables go to infinity, 
it cannot have new states coming out when ^/s increases. Spin structure is not a serious problem. Although 
Froissart only derived that bound for scalar particles [4^, we can imagine it is also applicable for 2 
2 quarks /gluons. Yet we can not derive it out, just for a lack of some systematical ways parallel to 
Mandelstam representation. 
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(KN) parametrization ^9] 



(2-2) 



A + B\ogs + C log^ s 



(2) 



tot 



and the Donnachie-Landshoff (DL) parametrization (Regge fits) j50] 



(2-2) 



Ds' + Es 



-s 



(3) 



tot 



We choose them instead of arbitrary parametrizations simply because (i) some unitarity 
constraints or physical backgrounds have already been considered in their derivations, 
which may make them more reasonable. The other reason is that (ii) the experimental 
result oi p + p cross section for all subprocesses a^^^ estimated from proton-air inelastic 
cross section a^^'^^^ [511 ESI [53l [Ml ESI [56l [571 [SBj up to ^/s < 10^"^ eV in cascade processes, 
has also been widely fitted by such kind of parametrizations 0. Because cr™ / crf^l^^ is an 
important parameter in latter calculations, it will be convenient for us to compute them 
from the same kind of parametrizations. 

The continuations of dominate subprocess g+g — >■ gl+gl by the KN and DL parametriza- 
tions, are shown in Fig. [2] and [31 However, some qualitative analyses show that, usually 
Eq. ([3|) can only construct mild concave functions in log-log diagram, just as crf^^^ in 
TeV regions. The reason is that when both D and > 0, the two terms in Eq. ([3|) 
are two asymptotes with a plus sign in between. Eq. ([3D can also construct convex func- 
tions in a log-log diagram when D ■ E < 0. That is, when the negative term becomes 
dominate, the curve will drop rapidly, which is not consistent with our purpose. Hence 
we simply fit cr{^/s) = Ds'^ in Fig. [HI with the physical meaning of retaining the term 
analogical to pomeron exchange but dropping out the term analogical to p, u, f and a 
exchanges. We notice that both the KN and DL parametrizations cannot do well to fit 
a^^{g + g gl + gl), which hints us that p + p (the KK states) are tremendously differ- 
ent from p + p — > (the (SMhadrons), hence appropriate parametrizations with reasonable 
physical background need to be explored. However, in this context, we use only our rough 
continuation for estimations below. 

2.3 Collision Producing 

The KK bosons and fermions can be produced hj p + p collisions for high energy protons 
accelerated in the Fermi mechanism [5^1 [SD], or a special lineage of it as the diffusive shock 
acceleration model [511 [521 [S31 [Ml [SSI [SSI [SZ] (see Ref. [SSl[Sn] for a review), with the cross 
section cr^f already been calculated in §2.11 and §2.21 We want to know how many of them 
can actually be produced in such environments. 

"^Extra subprocesses of p + p cross section from new physics beyond the SM, will hardly affect the 
applicability of old fitting (in the air shower experimental detectable region up to 10^* eV), because their 
branching radios are always neglectable. When the energy becomes higher, they may become more and 
more important, as shown in Fig. [2| and [3] However, we neglect this secondary correction in this context, 
simply because the large theoretical uncertainties. 
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Figure 2: Fitting of dominate subprocess g+g 91+ di by the KN parametrization. a^^/ is 
fitted in [58] also by the same parametrization, using the data of accelerator experiments, 
Akeno, Fly's Eye, Nikolaev and GSY [SH ESI ESI [56]. The data points in the top left 
corner are calculated from cxfotj"*^ in air shower experiments. The thin dashing lines are 
cross sections for subprocesses, as in Fig.[Tl and a™ is shown in a thick black line. Because 
we can not in fact do a beautiful global fitting of <7^^{g + g ^ gl + g^), as shown in the 
figure, we fit it emphasizing particularly on (i) a smoother derivative at the endpoint, 
or (ii) a preferable global fitting. We will always use the first strategy in the follow-up 
calculations in this context. 

To leave the pertinent computations for specific astrophysical sources to the later pub- 
lications, and to give a more universal applicable estimation, we discuss by imitating the 
same assumptions in deriving the Waxman-Bahcall (WB) bound as in Ref. [70l[7l]. Let np 
be the number density and Ep the energy of protons, the assumptions are as follows: (i) A 
percentage [fl of r] of the protons have an injection spectrum of dup/dEp = K ■ E'"' isotropi- 
cally distributed from their rest mass energy rrip = 938 MeV ~ 10^ eV to some cutoff energy 
-E'p,max in the source comoving frame, while a percentage of (1 — rj) of them stay at rest in 
that frame. In this paper, we will assume rj = 1 for simplicity. The generalization of our 
result to ?7 7^ 1 cases are straightforward. The power law spectrum in energy is a natural 
deduction of the Fermi mechanism ^9] or diffusive shock model [68] . However, for diffusive 
shock model, the existing derivations may need the assumption of collisionless or scatter- 

^The reason is that protons may need a threshold velocity v^nin to be actually accelerated. For example, 
in the diffuse shock model, the needed energy is a few times of mpWg/2 for protons to have larger gyroradii 
than the shock thickness [63| ■ 
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Figure 3: Sameline styles have been used as in Fig. [21 but otherwise the DL parametrization. 
Because of some technical properties of Eq. ([3]) which cannot make it really looking like 
a^^{g + g ^ gl + gl) (see the context for detail), in fact we only fit the cross section by 

ing only elastically with infinitely massive objects here and there, which is unfavorable for 
our p + p producing. Hence, a discussion of energy spectrum in a shock with influential 
inelastic scattering may be needed for our purpose, (ii) Sources are in the critical optical 
depth Tpp = 1, and the optical depth is independent of protonic energy. When deriving 
the WB bound as an upper limit, optically thin sources to p + p and p + 7 reactions have 
been assumed. Hence, the meaning of these assumptions is that the WB bound has been 
saturated. When estimating the upper bound of neutrinos, p + 7 interactions may be more 
effective than p + p interactions, thus they are given more attention. Producing of the 
KK particles by p + 7 reactions is also interesting (see also discussions at the beginning of 
[|2]); however, it is beyond the scope of this paper, by the lack of ready-made calculations 
about amplitude-squared and cross section. Here, we only consider p + p interacting cases. 
Detailed discussions about the competition of p + p and p + 7 reactions to produce both 
neutrinos and the KK particles are given in §3 however, a wiser way is to discuss them 
in association with specific astrophysical scenarios [721 [731 HH US]- Hidden sources with 
optical depth Tpp ^ 1 [76l [77] are very tempting for our purpose, because they can produce 
more KK particles by collisions. However, a big problem is how they can be accelerated 
and get away from the source successfully. We will discuss this situation together with 
some details of acceleration mechanism in §3.2^ comparing it with neutrino observations 
in [|5l Discussions linking to more realistic astrophysical environments are left to later 
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publications. 

If we assume that normal p + p interactions which dominate a^^^ do not change the 
number of protons[§, after one time of collision, the number of the KK protons in proportion 
to normal protons should be 

^ ^K-' f- f- E,',E. X (4) 

in which we choose ^/s = ^j2EiE2 to neglect the effect of rest mass and simply choose 
6 = 7i/2. The resultant Up^^j^/up depending on -Ep,max and a is shown in Fig. HI The cutoff 
energy -Ep,max in the source comoving frame can be calculated by comparing the size of 
the shock and the duration of an acceleration cycle, or acceleration gains and synchrotron 
losses of energy (see Eq. (Ill) and ([6]) for details). Sources of UHECRs should have -Ep,max > 
10^^ eV in the shock comoving frame, because of the reason that we have already seen a 
couple of UHECR events with energy > 3 x 10^° eV [7H1 EO] in the observer's frame. Even 
if neglecting the energy losses by the Greisen-Zatsepin-Kuzmin (GZK) effect [801 E] when 
propagating, and dividing a Lorentz factor F ~ 300 of the source comoving frame, we still 
need at least -Ep, max = 10^^ eV. We may prefer a = Sv^/^v^ — v^) for the nonrelativistic 
shock cases [SH], a = {3p^-2f3^Pl + P^) / {P^- p^) -2 38/9-2 ~ 2.2 for ultra-relativistic 
shock limit [HZ], where (P^) and {P^) the upstream and downstream velocities {P 
factors) in the shock frame, a < 2 by some statistical reasons of sources [82] , or a = 2 to 
be consistent with the assumption of the WB bound [TOllTT]. We see that the proportion 
^pkk/^p is not sensitive to either -Ep^max (for a large enough -Ep^max, which is easy to achieve 
for an actual astrophysical environment) or the parametrization strategies we choose, but 
sensitive to the spectral index a. 



3 Accelerating and Propagating of the KK Protons 

For the reason that the KK protons take charges, and may have a not- very- short lifetime 
(see discussions in ^ID before decaying to the LKP 71, they can also be accelerated by the 
same mechanism similar to normal protons. The acceleration process is never bothered by 
the KK cascade decay. After being accelerated, they may also propagate through the space 
for a considerable distance before decaying to 71. If the decay processes are forbidden or 
suppressed sufficiently by other reasons, they may propagate as far as banging into the 
earth. Hence in this section, we discuss accelerating and propagating properties of the KK 
protons, and also propagating property of decay intermediate and final states. 



3.1 The Fermi Mechanism and Diffusive Shock Model 

While discussing the accelerated properties and bounds of an extreme relativistic particle, 
some time scales are very important. We will discuss them in the framework of the orig- 

^It is reasonable because of baryon number conservation. Other final states are either unstable (like 
strange particles) or have small branching ratios (like the KK excited states). 
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Figure 4: The ratio ^^pkr/""-? after one time oi p + p collision for each proton in a E~" injec- 
tion spectrum, depends on -Ep^max and a. The results by the KN and DL parametrization 
are quite similar, thus we cannot distinguish them in this figure. 

inal Fermi acceleration system [591 EOl, that is, particles are randomly accelerated by the 
electromagnetic turbulence. However, because of the deep connection between the Fermi 
mechanism and diffusive shock model, some of the results are also suitable for the latter. 
We will calculate these scales in the source comoving frame (for the turbulence behind the 
shock, it is just the shock comoving frame). To transform to the observer's frame when 
the source itself is relativistic, both the size of the shock R and the energy of the particle 
E should multiply the Lorentz factor F, e.g., R R/T and E E/T. 

The time scales are as follows: (i) The mean escape time t-^ = R/c measures the 
time scale to get through the accelerating region, (ii) The Fermi acceleration time = 
rjRi^/cP"^ [83] measures how quickly the particle gains energy, where i?L = E/eB is the 
Larmor radius, B is the magnetic field strength, f] is the Alfven velocity, and ~ 8/3 — 
50/3 is a factor determined by the turbulent system, (iii) Synchrotron timescale tgy = 

measures the energy losses of synchrotron radiation, where ctt = 
87rrg/3 is the Thomson cross section, m = rrip or m = rriKK — is the mass for 

both normal or the KK protons respectively, and is the KK excited quark mass. 

Successful acceleration needs both > tA and tgj > t^, thus giving two limits of maximum 
acceleration energy 

Ei!.L - ^ (5) 
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and 




w?c^(3 I Que 



(6) 



rrie V V^tB 



Because of the same charges but different masses normal and the KK protons take, Emlx 

(2) 

cannot distinguish them but -Emax can. Hence, to find out the KK protons by the higher 
energy they can achieve, we need Emlx > Emlx or the source to be sufficiently large, i.e.. 



To emphasize that up to a constant, the constraint > is equivalent to the assumption 
R > Ri^, which is a universal requirement for accelerations correlated to magnetic fields. 

Notice that the formation of the cosmic-ray spectrum depends and only depends on (i) 
the average gain in energy per acceleration event, and (ii) the acceleration events a particle 
may suffer. For the original Fermi mechanism ^59j, the gain of energy is 6E = {v/c)'^E per 
collision in average (hence the energy should enhance exponentially), where v is the velocity 
of the reflecting obstacles (e.g., the electromagnetic turbulence). Hence, the spectrum index 
a = 1 + (c/f )^Ateia/Atineia dcpcuds ou the duration Ateia between elastic scattering (with 
infinite massive reflecting obstacles) and the duration Atineia to break down the energetic 
particle. For the diffusive shock model in non-relativistic shock wave cases [69] , the average 
momentum (hence energy for relativistic particles) gain for isotropic distributed particles is 
6p = 4:p{v^ — v^)/3vp, and the probability of escape per shock crossing is 4t>d/fp, where Vp is 
the velocity of accelerated particles. To generalize this argument to relativistic shock wave 
cases, isotropic distributions of particles are no longer applicable, thus an intuitionistic 
derivation is absent. However, we may assume that the qualitative phenomena are similar. 
Assuming that the total inelastic cross section for the KK protons is similar to that of 
protons, we have Atineia.KK ~ Atjnda.p, thus both (i) and (ii) are independent of particles' 
properties. The only difference between protons and the KK protons is the distinctness of 
their initial energy (rest mass). Hence, we have reasons to believe that the energy spectrum 
of the KK protons should shift to higher energy than protons by an amount of rnKx/'^^p- 

The resultant np-^^/rip for interzone of fixed observed E is shown in Fig. [5l which is 
independent of E for mxK < E < E^^^k- Notice that n^^^ is enhanced by a factor of 
(mKK/''^p)° by the spectrum shift, where a is the index of the spectrum. We choose 
a = 2.7 or 3.0 for the overall cosmic-ray spectrum observed below and above the "knee" 
energy -Eknee — 10^^ eV, and a = a for some special sources. For reasonable astrophysical 
sources, ^pkk/^p be not too small a number, to that we have opportunities to discover 
them by air shower identification. 

3.2 p + p Optical Depth Reexamined 

Taking no account of the details of especial astrophysical sources, the p + p optical depth 
Tpp (which is assumed to be 1 in §2.31 and subsequently quantitative estimations) can be 

^If the cross section for the KK protons is smaller than that of protons, wc will have the more energetic 
KK cosmic-rays to be observed. 




(7) 
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Figure 5: ^pkr/""-? fixed E when a = 2.7, 3.0 (for diffuse flux) and asource (for special 
astrophysical sources). It is not too small to be identifled. 

discussed by focusing on special acceleration mechanisms. 

In the case of the original Fermi mechanism [59l [60], Tpp > 1 can be achieved by a 
sufficiently long existent duration of acceleration sources. However, high energy protons can 
also be broken down by p+7 interaction or other processes, so several times of accelerations 
may be needed. 

In the case of diffusive shock model, things are a little more troublesome. As we have 
already indicated in §2.3[ the derivation of power law spectrum itself needs the assumption 
of coUisionless. In the mainstream tactics for particles to cross back and forth the shock 
front, the downstream particles are scattered by strong magnetic turbulence behind the 
shock, and the upstream particles are scattered by the self-excited Alfven waves [68]. The 
p + p collisional disintegration has to be competitive with particle missing in downstream 
scattering, which makes a constraint for its feasible parameter space (because we already 
know some restrictions from spectrum index, etc). However, we should not be too serious 
in that problem, because until now there is still no consistent computation in diffusive 
shock model to generate the kind of waves, scatterings and acceleration phenomena [57] . 
Scattering with analogously energetic particles (rather than with infinite massive objects) 
is also a good way to change direction but avoid loosing too much energy, especially in the 
upstream regions where self-excited Alfven waves themselves are in trouble. 
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3.3 The Modified GZK Cutoff 



For the case that (i) the KK protons have too long a hfetime to decay while propagating, 
when discussing the propagating of protons versus them, the most (or the only) important 
issue is the modification of the GZK cutoff [HOllHI]. Other issues while propagating, such 
as e^e~ pair production, cannot change tremendously the property of energy spectrum. 
The GZK cutoff, which is the threshold energy for protons to interact with the cosmic 
microwave background (CMB) photons to produce pions, as 



for m ^ m,r, is enhanced for the KK protons by mKK/''^p times. We have assumed 
that the quarks and gluons in pions are all the KK zero-mode (because the threshold 
energy to produce the KK excited pion is another rnKx/'^Tr times higher), and the cross 
section for 7cmb + Pkk — > Pkk(^kk) + 7r°(7r^) is still large enough to make a cutoff. The 
second assumption is irrelevant because the more observable KK excited cosmic-rays lend 
themselves to discovery; however, the number of particles is reduced tremendously for 
higher energies so we have little opportunities to really meet one. 

For the opposite case that (ii) the KK protons have already decayed to the LKP 71 (for 
other cases, including the decay intermediate states, the discussions are analogous; thus 
we will not discuss this further in this paper), they may also interact with CMB photons, 
as 7cMB + 7i X. However, we do not know what X really is, whether the cross section 
is sufficiently large to make a spectral cutoff, or what threshold energy these processes 
correspond to. Because the LKP 71 is an elementary particle (rather than the KK protons 
which include normal quarks/glouns), its cross sections may be much different from that 
of normal ones. If the process 7cmb + 7i — ^ 7i + vr° is sufficiently important, the threshold 
has an enhancement of rriKK/'^TT thus is much larger than the GZK cutoff energy, similar 
to Eq. dHD. 

4 Air Shower Identification 

If we can distinguish air shower events in which primary particles are protons/ions or 
the KK particles (the KK protons, the LKP 71 or other intermediate KK excited particles 
which bang into the earth), we are capable of giving a tighter bound on KK contamination. 
Otherwise, we can only constrain it when the KK particles dominate the mass composition. 

The mass composition of cosmic-rays is rudimentary, partly because cascade processes 
with the same initial condition are not identical with each other, partly because of our insuf- 
ficient knowledge of parton distribution, thus different cascade models (e.g. QGSJET [H] 
and SIBYLL [85]) cannot make consistent predictions with each other. An incomplete set 
of parameters to discriminate different compositions of cosmic-rays are as follows: (i) the 
elongation rate/shower maxima X^ax [Mj, when the energy of the particle can be decided 
separately by fluorescence technique, (ii) the magnitude of the fluctuation in depth of maxi- 
mum Xjnax [87], (iii) the number of particles nmax at Xmax, (iv) the speed of rise in nmax [88] . 
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(v) fraction of muons in the shower events, and some geometrically-base parameters like 

(vi) lateral distribution functions (LDFs), (vii) the thickness of the shower disk or (viii) 
the shower front curvature (see Ref. [89] for a review). Some of the diversities rise by the 
simulation results of cascade models [HH [HH]. Hence, to identify a KK excited cosmic-ray 
event for our purpose, the PDFs for protons take the KK charge and corresponding air 
shower simulations of (i) Pkk, (h) 7i or (iii) other intermediate KK excited particles are 
needed. However, they are beyond our scope in this paper. 

Naively, differences mainly rise for the reason that ions are made up from constituent 
nucleon of relatively lower energy (typically as (iii, iv) [HE] and the Zatsepin effect [90| 191]) 
cannot be used for reference, so a careful filtration of the derived parameters from above is 
needed. For some parameters, their values can only be determined statistically (such as (i) 
in Fly's Eye data processing [86]) at present, and are thus not applicable for our purpose. 

For the case of (i) the KK protons, similar total cross sections cxtot for p + p and p + Pkk 
are expected, because pkk has two/three normal valence quarks (and uncountable sea 
quarks) sharing energy, and atot {p + p) is not supposed to increase tremendously for larger 
center-of-mass energy ^/s (see the discussions for pkk production in high energy p + p 
collision in ^ for detail). Hence, the primary KK protons should also trigger a priori 
cascade processes. Differences may rise from the mass hierarchy between p and Pkk? ot 
equivalently from the energy transformation between normal and the KK quarks/gluons 
that should absolutely happen by color confinement, thus making the differential cross 
section da{p+p) / dQ and da{p+pKK) / dVL dissimilar. That may make the first few collisions 
much different; however, the phenomenology for the secondary showers of cascade particles 
like TT^, 7r°, p, p, and 7 are the same. Some similar work in Ref. [92] and [93], in which 
the cascade process of hadronized gluinos (analogous to our KK protons with the KK 
quarks/gluons in them) bang into the earth is simulated, indicate that the identification is 
possible; however, other researchers are suspicious of their result [9l]. It will be favorable 
for us, if the former result is correct, and can be naturalized to our KK proton case. 
Observable air shower events for the KK electrons are also expected, because it has turned 
on quantum electrodynamics (QED) interaction, thus the cross section should not be very 
small. However, the origin of the KK electrons is not discussed in this context. As a 
result, we can only distinguish primary protons/ions or the KK charged particles for their 
different cascade properties, and quantitative analyses are needed for the identifiability. 

For the case of identifying (ii) the LKP 71, things might be a little easier. We have 
already had trustworthy methods to identify photons from air shower data [951196] . because 
they interact with the geomagnetic field, thus starting the cascade process much before 
that of protons. Hence, if 71 can also interact with geomagnetic field, we can identify them 
easily. However, the quantitative calculations need detailed properties of 71 interactions. 

Exotic cosmic-ray events have already been regarded, for example in Ref. [97] . However, 
they may be particularly interested as the weakly interacting low energy ones, which have 
their first collisions inside our scientific equipments. For our purpose, it is more pertinent 
to find UHE exotic air-shower events. If they would happen, they may be explained by the 
UEDs scenario we mention in this paper (or other new physics models), and their absence 
can constrain the same theoretical models. 
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If we cannot identify the KK charged particles' cascade events from normal protonic/ionic 
ones, we can only constrain it when the KK particles dominate the mass composition. A 
cumbersome issue is that energy measurement is also related to the assumption of primary 
particles' composition. So that even if there have already been the super-energetic KK par- 
ticles (with energy largely exceeding the GZK threshold) banging into the earth, we might 
not discover them by our energy estimation. A wise way for energy measurements is to 
choose some kind of calorimetric measurements (like fluorescence light emissions [98| 199]) 
which are relatively model independent. Here in this paper, we assume that energy spec- 
trum can be determined without overall departure for protons/ions or the KK particles. 
Thus we are capable of discovering them when they dominate the mass composition, lead- 
ing to some identifiable phenomena such as a lack of the GZK cutoff, etc. 

As shown in Fig. [5l if assuming Tpp = 1 and a = 2, the resultant KK contamination 
npj^j^{E)/np{E) may be as large as 10~^ to 10~^. Hence if we can identify them individually 
from protonic/ionic ones, it is not difficult to make a discovery. However, if we can not in 
fact identify them from air shower data, a larger Xpp may be needed. 

5 Tpp Constraints from High Energy Neutrino Detec- 
tors 

The p + p optical depth Tpp is a crucial parameter for the applicability of our methodology; 
however, it should be discussed in relation to special astrophysical sources. Quantitative 
calculations are difficult because the models for UHE particle origins themselves, such as 
gamma-ray bursts (GRBs), active galactic nuclei (AGNs) or supernovae (SNe), are also 
full of dubious issues. 

However, an universal estimation can be made from the observation of neutrinos. The 
reason is that p+p collisions can produce the KK protons, but also have branching ratios to 
produce tt^, which decay to neutrinos mainly in modes tt^ ^ fi"^ + ^fii^fi) —>■ e"^ + i^ei^^e) + 
i^tj. + i^fj. B Thus if we have known (or have an upper limit of) neutrino flux in earth, we 
can use it to restrict optical depth Xpp as well as the KK cosmic-ray flux. Of course, the 
constraint can only be used as an upper limit, because we do not know what percentage 
of neutrinos are in fact produced by p + 7 interactions compared to p + p interactions, or 
whether the KK particles can be accelerated successfully, or whether they can leave the 
source as easily as neutrinos do. 

There are already some publications of p + p collisions of the Fermi accelerated protons 
in astrophysical environments [721 E31 EH ESI ESI EH HOI] - However, they are not of 
much use for our purpose, mainly because: (i) Their calculations are always oversimplifled, 
because p+p interactions are always unimportant, compared to p+7 interactions in speciflc 
astrophysical environments, (ii) Their discussions do not include subprocesses opposite to 

^Other processes such as — * + z^e(f'e) + t^fj, + i^fj. may also be important for our purpose, be- 
cause kaons loose less energy before decaying into neutrinos, by their relatively larger mass and shorter 
lifetime |100j . However, we will not discuss this further in this paper. 



17 



p + p 7T + X. (iii) Within some of their scenarios, one of the reacting protons does 
not suffer Fermi acceleration. In our scenario, it is efficient (and also favored) that both 
protons are accelerated, (iv) They make few attempts to calculate neutrino spectrum, or 
still assume ~ 0.25Et^ ~ 0.05Ep all the same. The latter simplification is reasonable 
in the p + 7 cases, but not legitimate for p + p cases, because usually Ep ^ nip ^ E^ 
makes Fp ~ F^r ~ {Ep + E^)/y/s, thus energy can be shared roughly by mass proportion. 
Although spectra of thermal p + p interaction have already been discussed in Ref . |102] ; 
however, a careful calculation of energy spectrum of nonthermal p + p collisions (especially 
in the case where both protons are accelerated by the Fermi mechanism) is absent. Briefly, 
only when both protons have roughly the same (direction and amplitude of the) momenta, 
the simplification E,^ ~ 0.05Ep can be used; however, the spectrum of neutrinos should 
be suppressed by rather than E'"" in that case. Hence, we may guess the energy 

spectrum of p + p neutrinos is steeper than E''^ in higher energies. Some numerical 
results for specific astrophysical environments show that the energy spectrum looks like 

^-(a+l) [73] ^jii 

give a more detailed universal calculation in §5.21 
5.1 p + p ^ TT + X Cross Sections 

The cross sections to produce vr^ mp + p^n + X collisions have already been calculated 
in |103lll04lll05[ll06] (see also |107[ll08j for recent developments). These parameterizations 
are mainly based on relatively low energy terrestrial experiments, and focus on Lorentz- 
invariant differential cross section (LIDCS) E{d^a/d^p) rather than total cross section 
(Ttot- Simple estimations show that they may behave well near the center-of-mass threshold 
energy ^/s > 2mp + m^ to produce pions, but increase too quickly to be consistent with the 
total cross section for all p + p subprocesses cxfolj^ [19], [50l EH] we use in this context. Here, 
to give a consistent but not too cumbersome estimation from an astrophysical (rather than 
collider physical) viewpoint, we use Badhwar parameterization |105j when cr'^^+cr'^ +0"^^' < 
'^tot' ~ 40 mb (thus ^/s < 2.98 GeV), and assume all af^^^ produce pions for larger ^/s with 
total cross section a^^^ for a fixed proportion as the same as when a'^^ +a'^ +a'^ = cr^ot'. 
The result is shown in Fig. El In the calculation, we assume (i) each p + p interaction 
produces single pion, and (ii) nix — 2mp [109] when evaluating the P||,max in Feynman 
scaling variable xp = p\\/p\\,raa.x- The applicability of both assumptions must be considered 
carefully for our purpose. 

For a more accurate estimation of p + p neutrino spectrum, the LIDCS E{d^a/d^p) 
(which depend on ^/s, the pion energy E* and the scattering angle 6* in the center-of-mass 
frame 0) should be used directly to evaluate F^ = Fcmr^(l + ■ v*/c^). However, it 
pours oil on the flames in the further complicated calculations. Hence, we integrate out 
E* and 6* by 

a-{V~s)=2n sinW/ dp^ x ^^== ■ E^ ( -—] , (9) 
Jo Jo y/pf + ml V« P\/ 

^In this papaer, we use superscript * to denote the center of mass frame, and normal characters for the 
source comoving frame. 
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and memorize p* thereby T* = E^/itlt, = y/p*2 _|_ m'^/m^ only by a weighting average 

{p:){V~s) = JpU---)/ /(•••) (10) 

where /(•••) = (t'^{\^). In addition, we neglect the effect of ■ v*/c^. The result is in 
the top left corner of Fig. El We see that for ^/s not too close to 2mp + m.^-, it is always true 
that y/s — 2mp — m^r ^ (p*). However, the result is not suitable for larger center-of-mass 
energy, and we lack a reasonable extension for ^/s > 2.98 GeV. 




2.0 2.5 3.0 3.5 

■/7 (GeV) 



Figure 6: Numerical result of total cross section a'^{^/s) within Badhwar parameter- 
ization [105] . The dark region shows that the parameterization is only suitable for 
^/s < 2.98 GeV, thus we force-fix the proportion a'^^ : o"'^ : a'^ for large ^/s. The 
weighting average {p*){^/s) is in the top left corner of this figure, with the same abscissa 
as cr'^(-yi). 



5.2 Neutrino Spectrum 

Hence, if we know the spectrum and optical depth of protons (within some specific astro- 
physical scenarios), the cross sections forp + p to produce vr^ — + t'e('^e) + + i^fi, 
we can calculate the spectrum of neutrinos produced hj p + p collisions. However, the 
calculation is not easily done. Even if decoupling the relationship between vectors and 
V*, and using only a weighting average momentum (p*) for some specific ^/s, the pion 
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distribution function should look as complicated as 




where ■\/s{Ei, E2,6) is the center-of-mass energy in the source comoving frame. To deal 
with the Delta function in an easier way, we fixed = 0.2 GeV and integrate 

out 6 visually. However, our simplifications may underestimate the harder part of the 
spectrum. To estimate the absolute value of flux to normalize our calculation, we use the 
same assumption as in deriving the WB bound [TOl [71] , that is, the cosmic-rays' energy 
production rate in the nearby universe is (-Ecr dhcji/ dEcn) z=o = 10^^ergMpc~^yr~^, or 
equally E^-^ dncn/ dEcB. = 5.05 x 10"^ GeV cm"^ s~^ sr~^ (or 3 times larger if thinking over 
redshift evolution). 

The resultant spectra are in Fig. [71 Our simplified estimation of vr and z/^ spectra 
are consistent with the numerical result E'^^^+i) in Ref. [73]. In our estimation, we have 
already considered the effect of neutrino oscillation (detectors are only sensitive i/^ flux). 
The WD bound can only be treated up to a constant, because if using the same treatment 
as ours, they should be much lower. The KK flux is estimated inaa = 2,a = 3 
environment (thus we use n-^^^^E) / n^i^E) ~ 10~^), other choices of test parameters can be 
read out directly from Fig. [5] We see that Xpp is not under strong constraints by neutrino 
observations. However, when we have more advanced neutrino detectors, they may give 
stronger constraints on Tpp hence the KK flux. 



5.3 Other Possible Constraints of Tpp besides Observations 

Another possible constraint of Tpp besides neutrino flux observation is GeV — TeV gamma- 
ray flux observations. The reason is that p + p collisions can not only produce t^"^ —>■ 
e'^ + i'e{i^e) + ^fi + i^fi, but it can also produce 7r° 2j. Hence, the EGRET, GLAST or Swift 
telescope can give it a constraint. Ref. [10]|, in which one energetic and one rested proton 
collide to produce a not-very-short-lived strongly interacting massive particles, considered 
both neutrino and gamma-ray flux constraints. We lack of discuss the latter one, because of 
the fact that both (i) to avoid this paper to be too long (the discussion is very similar), and 
(ii) to be different from energetic neutrinos which can only be produced by violent or high 
Fcm collisions, energetic gamma-rays can also be produced by the more mild synchrotron 
self Compton (SSC) processes |1T|. Hence the gamma-ray constraint is not as naive to deal 
with as the neutrino ones, so we leave it to later publications. 
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Figure 7: z/^ flux from p + p ^ tt + X collisions. The shadow regions denote cosmic-ray 
flux, with the KK ones for Tpp = 1 and Tpp = 10'* respectively. The detector information is 
cited from |llUj . Current neutrino detectors can only give constraints that Tpp < 10^. 

6 Discussion and Outlook 

Phenomenologies link ambitious theoretical physical models to reality, thus make physics 
a science 0. Terrestrial experiments are one of the ordinary methods to constraint new 
physics models; however, they have limited power because of our finite energy sources on 
earth. Cosmology can open an extraordinary window for new physics studies; however. Big 
Bang (and the extreme physical environment it has) happened only once in our universe, 
thus makes re-enactment impossible. Astrophysical constraints always have larger scope 
than terrestrial experiments (because they do not have the upper threshold of maximum 
achievable energy); however, meaningful scenarios presently known are always only avail- 
able for weakly interacting light particles. The motivation of this paper is to search for 
another way to construct new physics beyond the SM. 

In this paper, we construct an astrophysical scenario to (dis) confirm new physics for 
heavy particles beyond TeV energy scale. In our scenario, the KK protons are produced 
hj p + p collisions in Fermi accelerated environments, and they themselves are accelerated 
by the same environments. Because they may change the compositions and proterties of 
cosmic-ray events, air shower experiments can give a constraint to their properties. 

To know whether our scenario can give meaningful constraints to UEDs (and maybe 
other new physics models in later researches), we make some quantitative estimations. We 
first investigate whether enough KK excited states can be produced hj p + p collisions. 

^^At least in the philosophy of Karl Popper or previously. 
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We calculate the overall KK cross sections by precalculated Feynman rules and amplitude- 
squared, with also the CTEQ6.6m proton PDFs. Subprocesses g + g gn + 9n ^ + 1 ~^ 
may be most important in quark level. Because of color confinement, the KK excited 
quarks and gluons should form the KK excited protons. We then calculate what percentage 
of the KK protons can be produced by an isotropic and power law distributed proton 
spectrum. For the spectral index a = 2.0, Up^^^/np ~ 10^^° as an overall contamination 
after one time of p + p collision. However, when considering that the KK protons are also 
accelerated by the Fermi mechanism, the scene is much different. It is reasonable to believe 
that the energy spectrum of the KK protons should shift to higher energy than protons by 
an amount of wkk/'^^p, hence for some fixed energy E, the KK states should contaminate 
10^^ (for some special astrophysical sources with the spectral index a = 2.0) to 10~^ 
(for diffuse flux) of cosmic-ray events. So, if we can identify them from other cosmic-ray 
particles from air shower data, our method is capable of giving meaningful constraints. We 
notice that the GZK cutoff energy also shifts by a factor of mYj^/mp to higher energy, if it 
still exists. Hence, observations of cosmic-ray particles with energy much above the GZK 
cutoff, are given a reasonable explanation by the KK particles; however, the possibilities 
of this kind of observation is really small, even if the GZK cutoff does not exist. We 
also investigate the possibilities to identify the KK cosmic-ray events by air shower data. 
The investigation is still superflcial, because quantitative simulations are needed; however, 
the LKP 7i may be easy to identify, because they may interact with geomagnetic flelds 
just like normal photons, thus make the air shower tomography much different from that 
of protons. Finally, we calculate the possible constraints of the KK cosmic-ray flux from 
neutrino detectors; however, the constraints are very loose for current scientiflc equipments 
to affect our former estimations. 

It is appropriate to regard our calculations (outlined in this paper) as an "existence 
proof" for this kind of methodology. In fact, any charged particles beyond the SM, which 
are neither too light nor having a too short lifetime to suffer a Fermi acceleration, are 
suitable for our scenario. One immediate example is in UEDs, which we do not discuss 
in this paper because the extension is really straightforward. can cascade decay to 
Li or Vi [13]; however, because of their analogous masses, the lifetimes of Wf" should be 
much longer than in the Glashow-Weinberg-Salam theory of weak interactions. In 
fact, is also an intermediate state of our gi decay in this paper. Some lineage scenarios 
of the ADD model, which allow bulk bosons rather than bulk fermions [6l [71 [9], can also 
excite Kaluza-Klein which suffer our accelerating. However, a careful calculation 
of production rates and lifetimes is absent. There are also a lot of lineage scenarios of 
the RS model which allow the bulk SM flelds 11121 1113] . These models are more 
reasonable, because the RS model has an inherent orbifold conflguration (to obtain chiral 
fermions), and the bulk SM flelds can help us to understand some stiff physical problems 
like fermion mass hierarchy. Charged sparticles in supersymmetric (SUSY) models are also 

^^Of course, this scenario is also suitable for lighter particles. However, we can restrict the parameter 
space (of the endlessly emerging new models) tighter by other astrophysical/terrestrial methods. Hence, 
this scenario may specialize in new physics particles beyond TeV energy scale, especially just above the 
energy scale the best colliders can in touch. 
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good candidates for these kind of scenarios |114j . In order to explain the non-baryonic dark 
matter, the hghtest sypersymmetric particle (LSP) is preferred to be electrically neutral; 
however, it is not supposed to do so. Even if the LSP is really neutral, they can also be 
accelerated by the Fermi mechanism if the lifetimes of the charged ones decaying to them 
are not very short. Whereas different from the case of the LKP 71, if the LSP is gravitino 
or the lightest neutralino, they may hardly cause air shower processes because of their 
relatively small cross sections, even if they bang into the earth with tremendous energy. 
Gluino (which may exist as the form of gluino-containing hadron, compare to our KK 
proton) in split supersymmetry is also a good idea. There has already been one paper in 
Ref. PO], in which the gluinos are produced by astrophysical p + p collisions; however, (the 
astrophysical aspect of) our scenario has a lot of advantages than theirs, include: (i) One of 
the protons in their scenario stays at rest, hence the center-of-mass energy a/s = a/ f^pEp of 
p+p collisions is at least 10^^~^^ eV, which is no more than p+ (air) (for UHECRs to collide 
with the atmosphere hadrons) center-of-mass energy here in earth; however, because both 
of the protons in our scenario are Fermi accelerated, the center-of-mass energy of our p + p 
collisions should be at least 10^^ eV (we have already derived an overall Lorentz factor 
F ~ 300), which is impossibility for any other scenarios to achieve near earth. Despite 
the fact that we do not know whether our ideas of (renormalizable) quantum field theory 
or new physics nowadays are suitable for such a huge center-of-mass energy ^/s, we know 
that something should happen there, (ii) The maximum energy a gluino-containing hadron 
can achieve in their scenario is only 10^^'^ eV [92J; however, the maximum energy of our 
exotic cosmic-ray particles, can be even larger than the GZK cutoff (see §3.31 for a detailed 
discussion). Hence, because of the fact that the cosmic-ray spectrum itself has a large 
negative power law index of about —a ~ —2.7 to —3.0 (below or above the "knee" energy), 
for gluino-containing hadrons of energy ii^exo produced by protons with energy Ep, their 
content of cosmic-rays with definite E has an additional inhibitory factor of (i?exo/-£'p)", 
which has an order of magnitude of (lO^^-^/Etnee)^ '' x (-E^kneeZ-^czK)^ ~ 7.7 x 10"^^. If most 
of the UHE protons did not produce gluino-containing hadrons, the inhibition should be 
even stronger. Thus even if the gluino-containing hadrons are recorded by our scientific 
equipment (e.g., the cosmic-ray observatories), they are very difficult to be found out by 
such a lot of events with similar energy. However, because our charged exotic particles (KK 
protons in the context) are also accelerated by Fermi acceleration, their content in the UHE 
region should be as large as 10~^ to 10~^, hence not very hard to be identified, (iii) We 
do not really need the charged exotic particles to be longeval enough to suffer the travel 
from the source to earth; it is enough that their longevities are long enough to suffer an 
astrophysical Fermi acceleration. Hence, we do not have to worry about some adolescent 
(comparison with the "baby universe" era) cosmological bounds, such as the predictions 
of Big-bang nucleosynthesis, or cosmic microwave background. Another very interesting 
particle candidate for our scenario is the charged massive particles (CHAMPs) |115] . It 
is an ambitious dark matter constituent as yet (hence it is longeval to suffer a Fermi 
acceleration). The association of our astrophysical scenario drawing in this context and 
the CHAMPs, is a very interest issue; however, we leave it to the later publications. 
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